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ABSTRACT 


An acoustic energy reflection experiment was conducted in a labora- 
tory tank containing water and an artificial kaolin sediment. The 
reflection coefficient for the sediment was measured for a limited 
range of angles of incidence at one frequency and one set of measured 
sediment mass physical and viscoelastic properties. A mathematical 
reflection model was derived in terms of the known sediment properties. 
The measured and the theoretically calculated reflection coefficients 
are compared. The measured reflection coefficient and a probable angle 
of intromission showed poor agreement with theoretically predicted 
values. For the sediment parameters of the relatively soft clay used 
here, the reflection model is quite sensitive to values of sediment 
density and sound speed and is only mildly sensitive to values of 
dynamic rigidity and sound absorption coefficient near the angle of 


intromission. 
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I. INTRODUCTION 


A. BACKGROUND 

The development and use of bottom bounce sonar and long range 
passive detection systems, and the use of sonar in shallow water all 
require knowledge of the effect the ocean bottom has on sound propa- 
gation loss. This effect is most often characterized by the quantity 
leetom loss (BL) in decibels, which is related to the reflection co- 
G€fficient (R), the ratio of the reflected acoustic intensity to the 
incident acoustic intensity, by the equation BL = -10 10946 Ree the 
large variations in bottom loss throughout the world ocean bottoms, 
coupled with the passive or active sonar figure of merit sensitivity 
to small changes in this parameter, dictate a need for accurate 
measurement of bottom loss characteristics in each area of the oceans 
where bottom sediment parameters are different. 

A direct method for collecting such information would be an 
oceanographic survey, of massive proportions, to determine the acoustic 
reflection coefficient at various angles of incidence and various 
frequencies of interest in each distinct bottom region of the world's 
Seeans., On the other hand, surveys of other physical characteristics 
of the ocean bottoms have been going on for years, with the result 
that some core samples and raw and analyzed data exist for practically 
every Oceanic bottom region. Although much information is available 
in the literature concerning the relation of reflection and transmission 


of acoustic energy to the properties of the reflecting and transmitting 





media, a round table of exploration geophysicists concluded recently 
that their understanding of the "nature of seismic reflection is far 
from complete."* In anti-submarine and submarine warfare, where the 
exact nature of the reflection and transmission of acoustic energy 

at the ocean bottom must be well known, models for calculation of re- 
Election coefficients generally are limited to consideration of a plane, 
homogeneous, compressional wave, transmitted at various angles of in- 
cidence in water and reflected from and refracted in layered sediments 
having various densities and other physical properties, including 
acoustic energy attenuation. 

The energy loss a plane compressional wave undergoes upon reflection 
from the ocean's bottom is examined in a number of reports, papers and 
textbooks. Three of these references were chosen to describe the nature 
of the problem examined in this paper. White [1] was used as a basic 
text which describes the energy loss as a result of both compressional 
and shear waves being refracted in the sediment. Hamilton [2] uses the 
physical properties of cee compressional and shear wave speeds, 
along with related values of bulk modulus (K) and modulus of dynamic 
rigidity (ut) in the sediment to predict values of reflection coefficient 
and bottom loss. Hamilton treats the sediment as both an elastic and 
a viscoelastic medium. Bucker et all3]compares measured and theoretical- 
ly calculated values of bottom loss for multi-layered viscoelastic 


ocean sediments. 





* 

Aisetey ss iemeweand Allen, S. J., “What Direction Should We Set 
for Hydrocarbon Exploration?" Exploration Geophysics, Vol. 21, No. 3, 
p. 412-421, September 1973 





A number of the acoustic reflection experiments that have been 
conducted at sea, such as those referenced in [2] and [3], have shown 
varying degrees of agreement between theory and experiment. The dis- 
m@eeeement Usually is greatest near the critical angle or angle of in- 
tromission. The theoretical modeling of these experiments has been 
complicated by sediment layering and other non-uniformities and by 
@arficulty in obtaining realistic physical properties. A laboratory 
experiment involving a smooth-surfaced homogeneous sediment and an 
undisturbed water layer helps to eliminate some of these complications 


in comparing measured and theoretical reflection coefficients. 


B. DYNAMIC RIGIDITY 

A key element in improved acoustic reflection models has been the 
inclusion of the dynamic rigidity and dissipation properties of the 
sediment. Most researchers, such as Hamilton [2] and Bucker et al [3], 
used measured Stoneley wave speeds to calculate the real part of dynamic 
rigidity. Only estimates of the imaginary part were available. Also, 
Stoke's assumption that in a viscoelastic fluid the bulk viscosity is 
equal to zero was employed by these researchers to estimate a value of 
Pie imaginary part of the first Pamewconstamt (i): 

At the Naval Postgraduate School, the torsional oscillator rethod 
has been used to measure the real and imaginary parts of the complex 
dynamic rigidity. Cohen [4] described one set of these measurements 
for a saturated, water-kaolin artificial sediment, and Martinek [5], 


conducting research on the same type of sediment, measured values of 


absorption coefficient. 








Gee REELECTION COEFFICIENT 

An objective of this research was to examine the usefulness and 
Validity of the values of dynamic rigidity and absorption coefficient, 
along with other data on the physical properties of the sediment, in 
calculations with a viscoelastic reflection model, and to compare the 
calculated results with experimental reflection coefficients. The 
Meeeeclastic reflection model of Bucker et al [3] was modified to avoid 
the use of the assumption that the bulk viscosity is zero for a compres- 
Sional wave. Some preliminary experimental measurements were made of 
the reflection coefficient at one frequency and for one sediment condi- 


tion for comparison purposes. 
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Ii. THEORETICAL REFLECTION COEFFICIENT 


Ae DESCRIPTION OF THE MATHEMATICAL MODEL 

Figure 1 illustrates a plane wave reflection model which was 
assumed to be an approximation to the experiment. The symbols and 
notation defined in Figure 1 were drawn from Bucker et al [3]. Refer- 
ring to this figure, ¢ and Y¥ are potential functions defined in the 
following paragraph and x is positive to the right while z is positive 
downward. 

A summary of the development of water-sediment interface equations 
in [3] follows. The wave equation is given by the relationship 


U) ; ; ; 
Vo + (~)* @ = 0, and its assumed potential function solutions are: (1) 


>". = A‘, exp i(kx + aokz - wt), (2) 
o'7) = AY exp i(kx + a,kz - wt), (3) 
o"., = ANS exp i(kx - a kz - wt), (4) 
4 = BY, exp i(kx + bj kz - Wt) C5) 
i 
A" = O 1 — 3 
1 ; B 1 0 


In water only the real values of k are considered for 


(ies 0. 
k=-——— .. (6) 


C 
WwW 


The texm C is the measured speed of a compressional wave in water, 


A’, and AY5 are the measured amplitudes of the incident and reflected 


waves, respectively, AY and Bt are unknown quantities, 0. is the angle 


of incidence, and ory and wee both equal zero. 
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Figure 1. Theoretical Reflection Model. 
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For the water-sediment interface two boundary conditions may be 
@pplied to the potential function: 
1. the vertical particle velocities (w) are continuous across the 


interface; 





ap , av 
OZ : Ox (7) 


W = 


ume svertical stress—-tensor elementS are continuous across the 


interface; 


: 2 2 
~(-pw7> -2u a 


ax- OxdZ 


) (8) 


ae, 
II 


ZZ 


ZX 


Following Bucker [3], three interface equations are derived as follows 
from the boundary conditions and treatment of the sediment as a visco- 


elastic half-space; 





’ ’ a " = t t 
a, (A 5 A x a, A 1 + B 1 (9) 
a) 2 t " ens 2 t ' 
P. () (A 5 ae iE >) Ht (b, 1) (A ? +2b, (B 1)! (10) 
Dane ee 1) Brae. (11) 
1 ie l 
where: 
Pp 6” 
2 Ji 
a = ——— - ] (12) 
1 Mt Zu) 
ie °° 
2 il 
b, = ‘ l (13) 








= a = horizontal phase speed in water 
Sin @. 
a 
a, = cot 0 
2 
p, = density of sediment 
P, = density of water 


The Lamé constants A and py are complex, 





dy = My) - ih") (14) 


a9 vy! -— sai! 
u a ee ny (15) 


B. SEDIMENT PARAMETERS 

The above model permits a solution to the wave equation in terms of 
sediment parameters X, U, P and c, of which A is the only term which 
does not have independently determined values. Cohen [4] measured 


values of dynamic rigidity, w'. and uw". for a range of sediment densi- 
1 g 


JL 


ties. Martinek [5] measured the absorption coefficient, ane and the 


compressional wave speed in the sediment, c at a sediment density 


De 
comparable to Cohen's and at frequencies in the range of experimental 
interest. These values are summarized in Table I of Section IV. Addi- 
ional sound speed and density values were measured as part of this 
experiment and recorded in Table II of Section IV. 


In the absence of scattering from the water-sediment interface, as 


in theory is related to A". and uu". and c_. is related to A'. and u' 


i i al iL . 


With these known quantities, the next series of steps transposes the 


equations to eliminate A and ae as independent variables and solves 


for a reflection coefficient in terms of sediment parameters Wis we 


Pir Cy: along with Cu O. and 0. in the water layer. 
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Seen LCULATION OF REFLECTION COEFFICIENT 

In the following development, the symbols are Similar to those used 
femsorenerdt ([G] but with a reversed sign convention for the complex 
Seeomentiad, The Sign convention used here is the same used by Bucker 


Seal (3). Definition of terms is as follows: 


a 


k = k! - ik" = = complex propagation constant 
a= a' - ia" = complex wave speed and 
W = angular frequency 
4 Ly 
+ emma 
Soa 
ei = peVe 


Py 


where K, is the bulk modulus and p.v. indicates principal value. 


: 4 
+ it “fb ~—_, — ft ~ 
From Hamilton [2], K, 5 Hy ds +- 2u, Sou that 


A. + 2y ae 
Q= pov. ee ° (16) 


The compressional wave speed in the sediment, Cyr is given by 


i 


A) 


1 
= Re = ; hi 
[ 5 } i) 
The compressional wave absorption coefficient in the sediment 


is 


u 
= = Wim [ ee : (18) 


The complex compressional wave propagation constant is 


k = cue coer ee eo = eee 
Ja]? Ja]? cy De 








For physical reasons, a 7? QO to ensure decreaSing values of the 
Ce 
Peeential function»with positive increasing values of distance along 


@eempropagation path. By the same physical reasoning, c, > 0. 


1 
Peegebraic Steps to solve for reflection coefficient, R, in terms of 
Wee wee a os and frequency, f, and independent of and ae are 


shown in detail in Appendix A. These steps yield the following set 


of equations: 











De ? ? 
ee Z Z-1 
R= | x | oa | (A-10) (A-11) 
i: 
2 
A+ A" 2a ee (b Sails 
2 - Dl 1 
“at, 20, ~ | aera a areas as 
2 2 Pp. (b, +1)wW 1 
: | j 
ay =H + iW (A-12) 
Ys F 
a = + 71G (A-15) 
peu, 
F = —————-- - 1 (A-16) 
(is ye (UH yo 
iL 1 
2 
. P4¢ at 
— (A-17) 
diy ye ce "A 
1 ik 
Zeit) 32 
= —- = = -18 
anmelen 2 4 TE : (A-18) 
a q 
ect 
Es 
qe 
J = | H+ 1 (A-19) 
Clad «2 
ee 
Die 
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Equation (A-11) for R then was solved using standard WATFORG on 
the IBM 360 computer at the Naval Postgraduate School. Values of O. 
were chosen from zero to T/2 and the physical properties of the 
sediment were taken from Cohen {4]) and Martinek [5]. Values for BL 
and R were computed over the range of experimental error measured in 
the physical properties cited in these two references, in order to 


examine them for sensitivity to changes in these values. 


Did 








Mii ee Es SsUeE MENT OF REFLECTION COEFFICIENT 


A. SUMMARY OF THE EXPERIMENT 

The experiment was conducted by reflecting a pulsed beam of sound 
from a water-sediment interface at various angles of incidence. The 
incident and reflected wave amplitudes were detected with a directional 
receiver of design identical to that of the transmitter. A sound fre- 
Guency of 120 kHz was used. The received signals were amplified, fil- 
tered and signal processed to improve signal-to-noise ratio. The 
reflection coefficient was computed from the ratio of the squares of 
the voltage amplitude of the direct and sediment reflected pulses, 


uSing appropriate corrections for spherical divergence. 


eee tHe APPARATUS 

The artificial sediment, a commercially available pure kaolinite, 
was contained ina large, rectangular tank, 2.42 m long by 1.20 m wide 
by 1.00 m deep. The leveled, saturated sediment thickness was 60 cm. 
The water depth was 31.5 ene 

Identical transmitting and receiving transducers were set up in 
the water and mounted on a horizontal rail above the tank. Each trans- 
ducer could be rotated through 360° independently in both the vertical 
and horizontal axes and moved horizontally and vertically in the tank. 
The two electrostatic type transducers, built by research physicist 
Donald F. Spiel of the Naval Postgraduate School, have an active circu- 
lar face of 6.0 cm diameter, consisting of an aluminum coated, one mil 


thick plastic film mounted over an aluminum electrode which is set in 
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Sleeger plastic disc. in Order to keep the air gap between the plastic 
film and the electrode uniform, equalization of pressure between this gap 
and an air filled space behind the electrode was achieved by small holes 
drilled in the aluminum electrode. The inside of the transducer was 
partially evacuated. At 120 kHz the 3 dB beam width is about 20°. 

The basic geometry of the experiment 1S shown in Figure 2. Referring 
to this figure, both the transmitting and receiving transducers were 
rotated to a selected angle of incidence, 0. 

Based on the choice of centerline height of each transducer being 
23.5 cm above the sediment layer, the idealized ray path Yy was computed 
for each oe and the transducers were positioned along the horizontal 


rail so that Yy a Min and xX. = x The horizontal placement also was 


iL a4 
chosen to ensure that the same region of sediment was ensonified for each 
case. 

The range of angles accessible for this experiment is limited by 
the dimensions of the apparatus and the necessity to keep the sediment 
reflecting surface in the far field region of the directional trans- 
ducers. The extent of the far field region was determined experi- 
mentally by plotting the directly propagated sound field pressure as 
a Etection of range. 

The directivity of both transducers at 120 kHz provided a suffi- 
cient reduction in off-axis signals to eliminate interference from 
Side wall and water-air surface reflections. This was confirmed experi- 


mentally. At angles above 70° a small amount of interference between 


direct and sediment reflected pulses could be observed, the effects of 


which are discussed later. 





meaeer surface 
depth 


PO wea ect 
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Figure 2. Experiment Geometry. 
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Cc. INSTRUMENTATION 

Figure 3 illustrates the transmitting and receiving equipment and 
@eeewrtry in block diagram form. This configuration is similar to 
that used by Martinek [5]. The Hewlett Packard 204C oscillator fre- 
Suiemeyewas Set at 120 kHZ within + 3 HZ and was monitored with a fre- 
Meme y counter. The General Radio 1396-A tone burst generator was 
eo produce two-cycle tone burst pulses with an approximate 10 msec 
delay between pulses. This delay was sufficient to prevent undesired 
reflections from disturbing the measurements. The Hewlett Packard 
467A amplifier was used to boost the voltage by a factor of ten before 
meewsignal went to a 1:3 transformer for further amplification. The 
polarizing circuits for both the electrostatic transducers are shown 
ierigure 4. 

The received signal was amplified either 20 or 40 dB by a Hewlett 
Packard 465A amplifier and then passed through Krohn-Hite high and low 
pass filters set to a 20 kHz pass band centered at 120 kHz. A Prince- 
ton Applied Research (PAR) model 160 boxcar integrator was used for 
Signal processing. The time delay feature of the boxcar integrator 
was used to determine the times of pulse arrival. The enhanced signal 
from the arriving pulse was then fed to a Moseley X-Y Recorder. Fig- 
ure 5 which was traced from original records represents typical received 
wave forms after integration and recording. Wave form points selected 


to determine zero to peak voltages are noted on this figure. 


Pee DATA COLLECTION PROCESS AND ERROR ANALYSIS 
A number of angles of incidence were selected and appropriate hori- 


zontal separations were pre-calculated. At each setting of horizontal 


Ze 









Sine Tone 
Wave Burst 
Beencrator Generator 






Amplifier 
4074 












1:3 
7 0 mains Dual 
Beam Prans,tormer 
Cf 
Time Oscillo- 
aoe scope 
‘Cc 





Potat iz mg 
Caimcurt 


BoOxXxcar 
faecegrator 











x-Y 
Recorder 


printer 









Amplifier 
465/ 





Lolaioie, 11g 
Cieenrt 






Receiver 





Figure 3. Instrumentation. 
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Spacing, measurements were recorded of the received pulses for both a 
Girect path between source and receiver (0. =") and Lor the sediment 
reflected path. Additional gain was required for the latter measure- 
ment. The receiver then was moved its full horizontal travel to 
determine that no other pulses were being reflected from below the 
water-Sediment interface. 

Referring to Figure 5, zero to peak voltage was determined by 
averaging the maximum positive and negative peaks of each recorded 
pulse. In all cases these peaks occurred in the third cycle of the 
pulse. Voltage output anomalies in the tone burst generator, trans- 
formers, filters, amplifiers, boxcar integrator and the X-Y recorder 
were assumed to be the same for both direct path and reflected path 
voltages and therefore were cancelled by division in the reflection 
Seerricient ratio. Also, because of this ratio computation, trans- 
ducer sensitivities, sound pressure level calculations, etc., were 
not required. 

mae log-log plot of aL path average zero to peak amplitude ver- 
sus horizontal Separation waS constructed to determine the minimum 
horizontal separation at which the pressure varies linearly with the 
er se range (l/r). This established the inner limits of horizontal 
Separation which in turn set the minimum angle of incidence at 
0. = 50°. Since the reflected path was longer than the direct path, 
this separation insured that the ensonified patch of sediment was also 
in the far field. The dimensions of the tank limited the maximum 
angle of incidence to 0. = 716°. 

The arrival time of the sediment reflected pulse occurred later 


than the arrival time of the direct pulse. The measured arrival time 
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differences between the two pulses were compared with time differences 
calculated using the idealized geometry ray path difference and the 
known speed of sound. Compatibility between observed and calculated 
delays confirmed detection of the sediment reflected signal. 

The uncertainty in the angle of incidence (0. ) on each transducer 
as about + 1.0°. Part of this is due to the precision with which the 
scale could be read and part arises in the assembly process of aligning 
the transducer face to the zero mark on the angle indicator dial. It 
wees ctimated that horizontal separations were accurate to + 0.2 cm. 

The height of the transducer above the ensonified sediment patch was 
eeeurate to + 0.2 cm and the sediment surface within the ensonified 
patch was level and smooth within + 0.1 cm. With a sound wavelength of 
about 1.2 cm, this scale of roughness in bottom level 1s considered not 
to have a significant effect on observed signal amplitudes. The un- 
certainty in measuring voltage amplitude is estimated to be 2% of the 
Mrllavolt scale in Figure 5. 

At angles of incidence greater than 70°, small but perceptible 
off-axis energy from the direct path was detected. Since the 
difference between reflected and direct path lengths decreases with 
greater angle of incidence, the arrival time difference between the two 
pulses Saye so that the leading edge of the reflected pulse over- 
lapped the tail of the direct path ee pulse. This resulted in 


pesctimated + 5% error in reflected voltage amplitudes above 70°. 


EB. THEORY OF MEASUREMENT 


Referring to Figure 2, it was assumed that the sound pressure ampli- 


tude varied inversely as the square of the range from the source and 








that the absorption in the water was negligible. The assumption was 
also made that the sediment surface was sufficiently smooth and homo- 
geneous so that the reflected rays behaved as if they had originated 
at the mirror image position of the source. 

mae sreceived intensity I, at any range x along the direct path, was 
determined by orienting the beam axes of the transducers toward each 
other and measuring signal amplitude A'. Intensity I is proportional to 
eo. The reflected wave intensity I, which is proportional to the 
square of the measured reflected wave amplitude A', is a function of 
the reflection coefficient R, Tie and the spherical spreading along 
the total reflected path length y = Yy = ve 2yy0 assuming that the 
reflection is specular. I, is the axid)l intensicy at Unite distance 
from the sound projector. Stated in equation form, this relation is 
ae = (AN)? = RI, /y*. Solving for R and substituting for I, determined 


from the direct path measurements gives the following: 


2n"y, \ 2 


Oe 
Bottom loss is calculated by 


BL 


Pi 





IV. PRESENTATION OF DATA 


Tables I and II, which list some Se the mass physical properties 
of the sediment, are self-explanatory. Table II records the speed of 
sound and density in the sediment for three different depths near the 
surface. Sound speed ratio c,/c, was measured with a sound velocimeter 
and densities were measured using the sampling device described by 
Martinek {5]). The data show both density and velocity gradients in 
the upper part of the sediment. 

Table III lists a series of eight combinations of the five measured 
sediment parameters which are used in the theoretical model. These 
combinations are designed to determine the model's sensitivity to para- 
meter changes with each combination given a sequence number for refer- 
ence purposes. Sequence (1) is based on Table I dynamic rigidities at 
secaensity of 1.35 Ti ae along with sound speed and absorption coeffi- 
etent at a density of 1.38 mayen The maximum positive experimental 
error from Table I dynamic rigidities is the change inserted in Sequence 
(2). The effect of only a density change from 1.38 avin jeep dies Asie niV/ 
aa is determined from using Sequence (3). Together, Sequence (2) and 
(3) test the model for sensitivity to density changes. Sequence (4) 
Srotes to the Vable I values of dynamic rigidity for 1.28 anya density 
and in Sequence (5) the measured sound Suen for a density of 1.28 gm/ 
ae in Table III is inserted. Together, Sequence (4) and (5) determine 
the effect of only a sediment sound speed change. Also, Sequence (5) 
is based on the measured data which is assumed to be the closest approxi- 


mation to the sediment conditions during the reflection experiment. 
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Sequence (6) and (7) are designed to test the model for large arbi- 
Mmenry Changes in dynamic rigidity. The density of 1.24 Suen is an 
assumed value of sediment density at a point closer to the surface 
than that at which a suitable sample could be collected. Sequence (8) 
tests the model for a large arbitrary change in absorption coefficient 
relative to that in Sequence (5). 

maoles IVa and IVb ee the values of R and BL for the eight 
sequences computed from Equation (A-11) for a range of angles of inci- 
dence from 14° to 83°. From these two tables sensitivity of the 
theoretical model is determined. The values from Sequences (1) and (5) 
are plotted in Figure 6 for comparison. 

Table V tabulates the data from the reflection experiment. The 
Values of BL in this table are plotted in Figure 6 along with calcu- 


lated values from Sequences (1) and (5). 
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TABLE I 


Physical Properties of Saturated Water-Kaolin Sediment 


Cohen [4] Sediment Properties for 20°C 


Density Dynamic Rigidity 
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Martinek [5] Sediment Properties for 19.9°C 
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Measured Sound Speed and Densities 
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Not measured 
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Measured at 18.8°C 
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Corrected to 20°C, based on = 1480 x 107 cm/sec 
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Combinations of Sediment Parameters 
Used in Theoretical Computations 
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Angle of Incidence (degrees) 


Comparison of Results. 
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V. COMPARISON OF DATA 


A comparison of the bottom loss calculations from the reflection 
model with the measured results from the experiment reveals somewhat 
poor agreement between the two. The experimental limits on the angles 
of incidence prevented measuring tee en cone eee in the 
Miernity Of the theoretically predicted values of 40° to 43°. Figure 6 
shows a possible measured angle of intromission in the vicinity of 65°, 
but the lack of experimental results below 50° leaves open the possi- 
bility of even higher measured bottom losses at the lower angles. An 
additional comparison was made with a liquid-liquid model described in 
Kinsler and Frey [7, p. 145], using only the densities and sound speeds 
of Sequences (1) and (5), which predicted an angle of intromission in 
the range of 26° to 27°. The theoretical magnitudes of bottom loss 
near the angle of intromission are 10 to 20 @aB higher than the values 
measured at about 65°, 

tescight combinations of parameter changes reveal a number of 
facts regarding the theoretical model and its application to the experi- 
ment performed. A summary of comparisons among sequences 1S contained 
in Table VI on the following page. The high sensitivity to density and 
sound speed changes shown in this table agrees with what is expected 
from the basic liquid-liquid model [7]. The Table VI comparisons also 
reveal that dynamic rigidity and absorption coefficient affect bottom 
loss magnitude only near the angle of intromission. The range of para- 
meter changes selected moe the sensitivity tests had only a small effect 


on the range of theoretical angles of intromission, as shown in Table VI. 
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sequence 


Numbers 


1 and 2 


2 and 3 


4 and 5 


6 and 7 


eeand 8 





TABLE VI 


Comparison of Sequences 


Parameter 


Change 


Small change in 
Gaynamic rigidity 
Density 

Sound Speed 
Large change in 


dynamic rigidity 


Absorption 
Coefficient 


om. 


Samet vl ty Or 
Equation (A-11) 
Sensitive only at 


the angle of intro- 
mission 


Very sensitive 

Very sensitive 
Sensitive only near 
the angle of intro=- 
mission 

Low sensitivity only 


at-anagle Of iantromis= 
sion 








VI. CONCLUSIONS AND RECOMMENDATIONS 


The data collected for one sediment condition at one frequency for 
a limited range of angles of incidence must be considered as preliminary 
in nature. Three extensions of data collection using the developed 
experimental method are needed to complete the experiments. A variation 
in sediment conditions is needed to measure the effect on reflection 
coefficient from known changes in density, sound speed and dynamic 
rigidity. Apparatus now exists to measure all of these parameters. A 
Variation in frequency also is needed to measure the effect of changes 
in absorption coefficient, which has been assumed to be the only fre- 
quency dependent sediment property. Measurements over a wider range of 
angles of incidence are needed in order to determine experimentally the 
existence of an angle of intromission. In order to increase the range 
of angles of incidence changes in the apparatus will be required. As 
an example, a deeper and larger tank will permit measurements at smaller 
angles of incidence. 

The low sensitivity of R to values of dynamic rigidity could be due 
to Me vange of dynamic rigidities used for the material not being repre- 
sentative. Because of the fact that the reflection model is quite 
sensitive to sound speed and density, there exists the possibility that 
a part of the discrepancy between calculated and observed reflection 
coefficients is due to inadequate accuracy for measured sediment physi- 
cal properties. 

The sediment densities and velocities measured during this experi- 


ment and earlier by Martinek [5] confirm that gradients in properties 
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MemexiSt in the upper portion of the sediment. Although the gradyvents 
Mey Nave a Significant influence on the observed reflection, this was 
not accounted for in the theoretical model used here and time limita- 
tions did not permit the adaptation of another model. One such model 
which does include sound velocity gradients in the reflecting medium 
has been developed by Morris [8] based on the theory of layered media 
reflections described by Brekhovskikh [9]. .The gradient models must be 
included for experiments in which the sound reflection process involves 
greater penetration of sound into the sediment than was probably the 


case in the experiment reported here. 
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APPENDIX A 


Algebraic steps to solve for reflection coefficient in Section II 


ane set forth in this appendix. 
' 2 A/S) -j 
(Oo) = xX-aiy = (x + y ) / e” = re o (A-1) 


Note: @ not related to 0. 





= ‘ | — 
ah 4 + 2u 1 (A-2) 
vy =i _ 2\1 7 (A-3) 
x = ee 4- fae (A-4) 


fe ae ure = (Gy aeee O72 a cine / 2) 


Q = (x)? 


The positive roots of the following trigonometric identities were 


substituted: 
cos 6/2 = + Genres U2 Penney oan cos 8/2) 1/2 
a= cr) 7* (rey at 7? - se a 1°77} 
1 2 Ve 


v2 a2 


xa Warne —- 1(r-x) ] 
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; — ; ae ae 
Rationalizing the inverse of the above to obtain 5 Yields the 


following: 


y oF 


wage of ieee 


r V2 


Q Ire 
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Substituting the real and imaginary parts into Equations (17) 


aman{i18) yields the following: 








72 
Re [= a (r+x) 
aye 
V2 
Tn [= ee) 
r V2 
c) — eee (A-5) 
vr ae 
_ V2 mE /r-x (A-6) 
p ‘2 


Algebraic manipulation of c, and a then was used to eliminate 
p 


1 


yr and state all relations in terms of x and y. 
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Manipulating further as follows: 
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2 ZO 2) ies) 
pale i Oe ae 
271f — oye+x 
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LD eee 
Fe cc a. 
7 ee ia, 
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ey ere 
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and since as and Cy are greater than zero, the positive root is chosen. 


Combining x/xr and y/r to obtain y/x yields the following: 


ap 


Sg 2 arta tte ta (A-8) 
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Cad 
_ (2 P)? 
t (TE 





; Kee ee: Zon , 
Substrtuting the definitions of x, y and © in Equations (A-2), 
(A-3) and (A-4) shows that. u',, es AY and An are related only to 


the groupings of known physical properties of Equations (A-7) and (A-8). 
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Solving the interface equations by adding Equation (9) and (11) 


and substituting the result into Equation (10) gives: 


2 2 | 
(2 1) Cooma): 





op, (=)° (A', + a",) = yp oper 
rae WS 2 1 2a, cl ose 
0? er 
= —— 7 
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pe | 2 en ra (cs 
1 = it at 
Uwe omits | i 1) 2a, 


Letting Z equal the right hand side of Equation (A-9) and defining 


intensity reflection coefficient 


A", 2 
R = nn ako} 
2 
ne + ye it Se ye TE 
2 2 Z 1+ SR 
such that a Os L+ VR _ ?, 
' = " = " ' S 
ah eal el /R 
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7-1 
then R = | rT | (A-11) 
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All quantities on the right hand side of Equation (A-11) are physi- 


Cal properties of the sediment or water and Z is a complex number. 


Formulating ay 


H and J are defined 


F and G are defined 


and b 
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P, { 


San 6.” 
si 


2 
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At this point the two Lamé@ constants A'. and Ms are contained only 


J 


in the terms H and J. Equations (A-13) and (A-14) are then restated in 


terms of x and y as defined earlier in Equations (A-2) and (A-3). 
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Substituting from Equations (A-7) and (A-8) then yields the final 


expressions for H and J. 
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